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H , Abstract. Enhancement of strangeness production has since long been proposed as 

a promising signal of quark-gluon plasma production. A convenient indicator for it 
is the Wroblewski parameter which has been shown to be about a factor two higher 
i in heavy ion collisions. Using a method proposed by us earlier, we obtained lattice 

QCD results for the Wroblewski parameter from our simulations of QCD with two 
light quarks both below and above the chiral transition. Our first principles based and 
^ ' parameter free result compare well with the A-A data. 
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1. Introduction 

Amongst the many proposed signatures of quark-giuon plasma in relativistic heavy ion 
collisions, enhancement of strangeness production JT] is a promising signal. Proposed 
in the very early stages of this field, it was based on simple model considerations 
like many other signals. Exploiting the fact, that the strange quark mass is smaller 
than the expected transition temperature whereas the mass of the lowest strange 
hadron is significantly larger, it was argued that the production rate for strangeness 
in the QGP phase, crQGp{ss) is greater than that in the hadron gas phase, o"i^G(ss). 
While this energy threshold argument for strangeness production in the two phases 
is qualitatively appealing, one has to face quantitative questions of details for any 
meaningful comparison with the data. Since the temperature of the plasma produced in 
RHIC, or even LHC, may not be sufficiently high for perturbative QCD to be applicable, 
estimation of strangeness production by lowest order processes like gg —>■ ss could be 
misleading. Indeed, it is now well-known that even for charm production, the next order 
correction to gg — > cc is equally large. 

A variety of aspects of the strangeness enhancement have been studied and many 
different variations have been proposed. One very useful way of looking for strangeness 
enhancement is the Wroblewski parameter [2j. Defined as the ratio of newly created 
strange quarks to light quarks. 



the Wroblewski parameter has been estimated for many processes using a hadron gas 
fireball model jH]. An interesting finding from these analyses is that As is around 0.2 in 
most processes, including proton-proton scattering, but is about a factor of two higher 
in heavy ion collisions. An obvious question one can ask is whether this rise by a factor 
of two can be attributed to the strangeness enhancement due to quark gluon plasma 
and if yes, whether this can be quantitatively demonstrated. We show below how quark 
number susceptibilities, obtained from simulations of lattice QCD, may be useful in 
answering these questions. 

2. As from Quark Number Susceptibilities 

Quark number susceptibilities (QNS) can be calculated from first principles using the 
lattice formulation. Assuming three fiavours, u, d, and s quarks, and denoting by /i/ 
the corresponding chemical potentials, the QCD partition function is 



f=u,d,s 

Note that the quark mass and the corresponding chemical potential enter only through 
the Dirac matrix M for each fiavour. Defining /io = + A^d + /^s and fis = fJ^u — f^d, the 
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baryon and isospin densities and the corresponding susceptibilities can be obtained as: 

_Td\nZ _Td'^\nZ 
V dfii ' '^^'^ V dnidnj 

QNS in (j2)) are crucial for many quark-gluon plasma signatures which are based on 
fluctuations in globally conserved quantities such as baryon number or electric charge. 
Theoretically, they serve as an important independent check on the methods and/or 
models which aim to explain the large deviations of the lattice results for pressure 
P(/x=0) from the corresponding perturbative expansion. Here we will be concerned 
with the Wroblewski parameter which we ^ have argued can be estimated from the 
quark number susceptibilities: 

A. = ■ (4) 

Xu + Xd 

In order to use to obtain an estimate for comparison with experiments, one needs 
to compute the corresponding quark number susceptibilities on the lattice flrst and then 
take the continuum limit. All susceptibilities can be written as traces of products of 
and various derivatives of M with respect to /x. With m„ = diagonal Xm's can 
be written as 
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(5) 
(6) 
(7) 



Here O2 = Tr M-^M^ - Tr M'^M^M-'M^, and Oii(m„) = (Tr M'^M'^f. The 
traces are estimated by a stochastic method: Tr A = Ylf=i RiARif^Ny, and (Tr A)"^ = 
2Si>j=i(Tr v4)j(Tr A)j/L{L — 1), where Ri is a complex vector from a set of A^^,, 
subdivided further in L independent sets. 
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Figure 1. Comparison of quenched and full QCD (left) results and typical continuum 
extrapolation results (right). 

Left panel of Figure 1 displays results [S] for the susceptibilities as a function of 
temperature in units of Tc, where is the transition temperature. Normalized to 
the corresponding ideal gas results on the same lattice, i.e, the infinite temperature 
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limit, results for QCD with two light dynamical quarks of mass 0.1 are shown 
as points whereas the continuous curves correspond to the results in the quenched 
approximation. The valence quark mass m^, appearing in ((SI)-©? is shown in the 
figure in units of T^. Note that in these two cases differ by a factor of 1.6 but the 
results for the corresponding susceptibilities as a function of T/Tc differ by a few per 
cent only. Encouraged by this, we investigated the continuum limit for the quenched 
case by increasing the temporal lattice size from 4 to 14, as shown in the right panel of 
Figure 1 for T = 2 Tc. The continuum results for QNS thus obtained in the quenched 
approximation are exhibited in Figure 2 for small m^. 
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Figure 2. Quenched QCD results in continuum limit (left) and comparison of the 
corresponding with RHIC and SPS experiments (right). 

The strange quark susceptibility was obtained from the same simulations by simply 
choosing m^/Tc = 1 (in both full and quenched QCD in view of Figure 1). Using (jlj), 
Xs{T) can then be easily obtained, and extrapolated to Tc by employing simple ansatze. 
The resultant \s(Tc) in quenched QCD is shown in the right panel of Figure 2 along with 
the results obtained from the analysis of the RHIC and SPS data in the fireball modeljS]. 
The systematic error coming from extrapolation is not shown but is of approximately 
the same size as the shown statistical error. The agreement of the lattice results with 
those from RHIC and SPS is indeed very impressive. 
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Figure 3. as a ftmction of temperature for full and quenched QCD. 

The nice agreement needs to be treated cautiously, however, in view of the various 
approximations made. Let me list them in order of severity. 
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• The result is based on quenched QCD simulations and extrapolation to Tc. As seen 
from Figure 1, the quark number susceptibilities, and hence As(T), are expected 
to change by only a few per cent. Since the nature of the phase transition does 
depend strongly on the number of dynamical quarks, a direct computation near 
for full QCD is desirable. We are currently making such a computation and have 
some preliminary results for full QCD with two light dynamical quarks for lattices 
with four sites in temporal direction. These are shown in Figure 3 along with the 
continuum quenched results for As(T) and the band for experimental results. While 
the emerging trend is encouraging, further exploration with varying strange quark 
mass, temporal lattice size (to obtain continuum results) and spatial volume is still 
necessary. 

• The experiments at RHIC and SPS have nonzero albeit small /i whereas the above 
result used /i = 0. Based on both lattice QCD and fireball model considerations. 
As is expected to change very slowly for small /x. This can, and should, be checked 
by direct simulations. 

• Lattice simulations yield real quark number susceptibility whereas for particle 
production its imaginary counterpart is needed. Assuming that the characteristic 
time scale of plasma are far from the energy scales of strange or light quark 
production, one can relate [6] the two to justify the use of lattice results in obtaining 
As. Observation of spikes in photon production may falsify this assumption. 

3. Summary 

Quark number susceptibilities which can be obtained from first principles using lattice 
QCD contribute substantially to the physics of RHIC signals. In particular, the 
continuum limit of Xu and Xs, obtained in quenched QCD, leads to As(r). Its 
extrapolation to appears to be in good agreement with results from SPS and RHIC. 
First full QCD results near Tc confirm this as well, although many technical issues, e.g, 
finite lattice cut-off or strange quark mass, need to be sorted out still. 
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